Global pointwise decay estimates for 
defocusing radial nonlinear wave equations 



Roger Bieli and Nikodem Szpak 



o 
o 

(N 



Max-Planck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, Golm, Germany 

(Dated: March 4, 2009) 

We prove global pointwise decay estimates for a class of defocusing semilinear wave equa- 
tions in n — 3 dimensions restricted to spherical symmetry. The technique is based on a 
conformal transformation and a suitable choice of the mapping adjusted to the nonlinearity. 
As a result we obtain a pointwise bound on the solutions for arbitrarily large Cauchy data, 
provided the solutions exist globally. The decay rates are identical with those for small data 
and hence seem to be optimal. A generalization beyond the spherical symmetry is suggested. 
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We consider a class of nonlinear wave equations 



(1) 



in n = 3 spatial dimensions where the nonlinearity f(4>) is of a defocusing type. It means that 
the nonlinear term has a repulsive action on the waves and focusing of waves is suppressed by an 
energy condition. Essential is the sign of the nonlinear term which is chosen such that the following 
conserved energy is positive definite 
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l\dt<P\ 2 + l\V4>\ 2 + F(4>) ) d A x 
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with / = —F'. Equation (pQ) has been intensively studied in the literature over a few decades, 
in particular in the case of a pure power nonlinearity f(4>) = — |</>| p_1 </>. Let us collect the most 
important results defining the context for thi s work . The global existence of C 2 solutions to the 
Cauchy problem has bee n shown by Jorgens Jor6l| in the e nergy subcritical case I < p < 5 and 
later by Grillakis Gri90l ] in the critical case p = 5 (see also Str88] for the spherically symmetric 
p = 5 case) while not much is known about the supercri tical ca se p > 5. For more references related 
to global existence we recommend the bo ok of S ogge |Sog95 |. Uniform boundedness of solutions 
for 2 < p < 5 has been proved b y Pech er Pec74f | . Uniform d ecay l/ t 1 ~ t and scattering have been 
proved for 3 < p < 5 by Strauss Str68| . Baho uri and Shatah |bS98[ | have shown that finite energy 
solutions decay to zero for p = 5 and H idano HidOll ] has shown scattering and decay to zero for 
2.5 < p < 3. Ginibre and Velo GV89I ] have shown scattering in the energy space for 1 < p < 5 
in various dimensions, but in n = 3 only for small data. Scattering results imply, in some sense, 
that the solutions behave asymptotically like solutio ns of the linear equation. (For more scattering 
results we refer to the monograph of Strauss Str89l ] .) 

Here, we go further and study the pointwis e behav i or of solution s for 3 < p < 5. Such results 
in dimension n = 3 exist only for small data |Asa8fl IST97L ISzpOSl ISBCR07| and are based on 
perturbation techniques which cannot be generalized to large data. We extend the tec hnique of 
conformal com pactification developed by Choquet-Bruhat, Christodoulou and others in CBPS83L 
Chr86l . lBSZ90j . The novelty is that we use a suitably chosen conformal transformation, adjusted to 
the nonlinearity. In this setting, we first show uniform boundedness of the transformed solutions in 
the precompact region and then use the inverse transformation to get an extra decaying (conformal) 
factor and thus a pointwise decay estimate for the original solution. We claim that this estimate is 
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optimal (for generic initial data) a s it is iden tical with the one obtained for small data wher e it has 
been shown to be optimal Szp08l . ISBCR07I ]. Numerical simulations done by Bizoh et al Biz08] 
also support this picture. 



The main result and the method 



In the following we study the Cauchy problem for the radial semilinear wave equation 







A<(> 



(3) 



with the initial data of compact support r £ [0, a p [, where a p > will be specified later 1 , and 
chosen at time t = 1 



4>(l,r) = 4> , d t 4>(l,r) 

and prove that the solutions satisfy for all t > 1 and < r < t 

C 

\4>(t,r)\ < 7 T7 r k 

' n ~ (l + t + r)(l + t-r)P- 2 
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with some constant C depending only on 4>q, 0i and p. 

Briefly, our method is based on a conformal mapping (u,v) — > (u,v) defined in the double-null 
coordinates u = t + r,v = t — r, 



1 



1 



u :- 



v :- 
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The conformal factor ft := r - multiplies the transformed solution cp(u,v) := $7 1 4>(u, v) which 
satisfies the transformed wave equation 



□0 + 
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u — v 



(_^)l/(p-2) _ (_£)l/(p-2) 
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in a precompact region of spacetime. There, we are able to show boundedness of some pseudo- 
energy flux what we further use to show the uniform boundedness of \(f>(u, v)\ < C. Finally, the 
inverse transformation provides the desired estimate 



( U ,v)\ <c -n = c ■ 
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u P-2 _ v p-2 
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u — V 



u ■ v p 



The power p — 2 in the mapping © cannot be increased, what would potentially lead to a stronger 
pointwise decay of the solutions, because then the factor multiplying the nonlinearity in ([?D becomes 
singular and our boundedness theorems cannot be applied. 

The estimate ([5]) is also optimal i n the sen se of compatibility with the small data case. For 
small initial data it has been shown in SBCR07I | that the solutions behave asymptotically, for large 
t and fixed r, like 



4>{t,r) 



tP 



— x + 0(t~P) 



1 Symmetries of the equation Q allow for mapping any compactly supported initial data onto the interval r £ [0, a p [ 
at t = 1. 
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where the constant C can be expressed explicitly via (fto and (ft\. Hence, the estimate ([5]) applied 
to the small data case gives the optimal decay rate. 

The assumption of spherical symmetry is essential for the method to yield optimal decay for p > 
3. On MxR 3 the obvious analogue of the map (u, v) — > (u,v) is no longer conformal, unless p = 3, 
of course. The reason is that radial and angular parts of the wave operator transform diversely 
and the resulting equation is no longer semilinear. Applying the transformation corresponding to 
p = 3, which is conformal, to an equation with nonlinearity of power p > 3 still gives a global 
l/(m;)-decay result, being however not optimal and independent of the actual power p. 



II. CONFORMAL TRANSFORMATION 



In the case p = 3, the radial wave equation ([3]), having the explicit form 

2 



df<j> - d;<t> - -d r (ft + 0=0, 
r 



is not only invariant under time translation, (ft(t,r) i— > (ft(t + a,r), scaling, (ft(t,r) i— ► X(ft(Xt, Xr), and 
reflection, (ft i— > —(f>, but also under a conformal inversion which is given by 

1 / — t r \ 

Ht,r)^ ^—^ <t>[j2—^I '{T^i )■ ( 8 ) 

This inversion maps solutions on the forward light-cone K + := {(t,r) | < r < i] to solutions on 
the backward light-cone K~ := {(t, r) \ < r < — t} of the origin and vice versa. Thus, establishing 
boundedness of a solution on K~ towards the future immediately implies pointwise decay estimates 
for the transformed solution on K + . In order for this to yield optimal decay in the case of more 
general nonlinearities it is necessary to also consider more general conformal transformations which 
is what will be discussed in the following. 

On (i, r) G R x R + consider null coordinates u := t + r and v := t — r. Then u and v are positive 
on K + and negative on K~ . For any p > 2 define, with respect to the (u, t;)-coordinate system, 
the map 



<5>:K + ^K~, (u,v) 
It is analytic with analytic inverse 



1 1 



yp-2 vP - 



$ 1 : K -^K + , (u,v) i-> ((-u) p~ 2 ,(-v) p-2 



Furthermore, there exists a positive analytic function £1 > on K + with the property that 

rQ = r o $ (9) 

holds on K + . If 77 = (it 2 — dr 2 denotes the Minkowski metric on R x R + then its pullback by $ 
satisfies 

which shows that mapping by $ is conformal. In the case p = 3 the map (ft >— > Q5>*^> = • 
o corresponds to the conformal inversion (jSj) since then = l/(t 2 — r 2 ). This conformal 
transformation is illustrated in Figure [TJ Note that with 

1 , r 

Op := 1 — G ]0,1[ 

2~ 
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FIG. 1: The function $ maps the forward light-cone K + to the backward light-cone K~ conformally and 
bianalytically. An interval / of the line {t — 1}, containing the support of the original Cauchy data, is 
mapped to a curve H. The future of the interval / is mapped into the future of the curve H. Initial data 
on H is evolved in the region S to yield new initial data on the interval J of the line {t = —1}. 



the compact interval I contained in {1} x [0, a p [ implies that J is compactly contained in {—1} x 
[0, 1[, where J is the intersection of the future of the curve H := <£(/) with the line {t = — 1}. 
Moreover, for p > 3 one has a p < 1/2, 

Suppose that h € C 2 (K~) is a twice continuously differentiable function on K~ . Then Q<&*h is 
a function on K + of class C 2 for which 

D(m*h) = (nn)$*h + n ^'^L^ mh > 

From relation ([9]) it follows that DO = such that 

n(m*h) = (p - 2) 2 (uv)- ( - p - 1) n$*nh. (10) 

Hence, if p > 2 and cfi E C 2 (K + ) is a classical solution of the radial wave equation 

D(/) + (Pl^- 1 = (11) 
on its conformal transformation ip := <E>*(Q -1 (/>) € C 2 (K~) satisfies 

+ 7^72)2 [**(n««) p_1 ]V'^r 1 = 0, (12) 

where denotes the push- forward by the diffeomorphism $. On the other hand, it follows directly 
from pop that if ip £ C 2 (K~) is a classical solution of the transformed equation (|12|) on ET~, the 
function (f) := Q$>*ip € C 2 (K + ) solves the original equation (jlip on ivT + . 
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III. EVOLUTION OF THE PSEUDO-ENERGY 



As mentioned introductorily, for equations of the form (pQ) there is a conserved energy ([2]) that 
in particular controls the spatial L p+1 -norm of the solution. One can no longer expect to find such 
a quantity for the conformally transformed equation (|12p or, more generally, for an equation 

IP— l 



□V + c(i, r)^|-0| p " = (13) 

with a non-negative function c of class C 2 . However, an analogue of the energy expression, although 
not conserved exactly, will turn out to be sufficient to prove boundedness of ip on the relevant region 
of K~ , provided the function c is monotonically decreasing in t there. This pseudo-energy will now 
be defined. 

Let ijj £ C 2 (K~) be a classical solution of (PT3l) on K~ for p > 2. Then the vector field £ given 

by 

l 



£ 



i(3«*) 2 + i(a^ v ( 



dt - r 2 (d t ipd r ijj)d r 



is continuously differentiable and 



div £ = r 1 



p+ 1 



□V> + cV#l p_1 

r 2 (d<c)|Vf +1 



p+ 1 



r 2 (d t c)W\ 



P+i 



holds on K~ . Assume furthermore that dtc < is non-positive such that the same is true for div£. 
This implies, recalling Figure [TJ together with the assumption that the support / of the initial data 
is compactly contained in {1} x [0, a p [, that 



0> / (div £)dt A dr = / i £ (dtAdr) = 

JS JdS JJUH 



{igdt)dr — (isdr)dt 



.1 
V 



I W)2 + I( 8ri/ , )2 + _l_, 



IP+i 



dr 



2 Ji 
p-2 



igdt) o $(1 + r, 1 — r) • 
j (iedr) o $(l + r, 1 — r) 



1 



1 



(1 - r)P- x 
1 



(1 +r)P~ 1 
1 



(l_ r )p-l (i + r )p-i 



dr 
(ir, 



where is denotes the interior multiplication with the vector field £. The last two integrals over the 
interval I depend only on the initial data (j)Q and 4>\ on / and are certainly finite. Thus 

1 



Eh 



o ■- 



c\il>\ 



p+1 



dr<C 



(14) 



for a constant C. 

The pseudo-energy Eq now controls light-cone integrals of the quantity c|^| p+1 according to the 
following Proposition [H a fact that will be essential for proving boundedness of ip on the region 
K := K~ n {—1 <t< 0}. 

Proposition 1. Let tp G C 2 {K) be a solution of (I13p with initial data supported on J satisfying 
the estimate (fT4"j) . If dtc < on K then 

/toWs /or any — 1 < to < 0. 
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Proof. It is useful to consider the functions 



E to (t) := f° (i £ dt)(t,r)dr 
Jo 

for any — 1 < t < to < 0. Then, for a fixed — 1 < to < 0, Et is continuously differentiable and 

E' t0 {t) = -{i £ dt){t,t -t)- [ [d r (i E dr)](t,r)dr + /' (div £)(t,r)dr 

Jo Jo 



<t -tf 



2 p+l 



t -t 

(t,t -t) + I (div £)(t,r)dr 
o 



holds for all —1 < t < to- Hence, it follows that 



/to 
is (dt + dr) (s,to — s)ds 



to 



(to - sf 



2 p+l 

to pto—t 



IP+1 



(s, to — s)ds 



E' tQ (s)ds+ / / (div £)(s,r)drds 
l J-i Jo 

E t0 {-1) + f° f ° (div S)(s,r)drds <E , 



since Ef (to) = 0, Et (—1) < Eq and div£ < 0. But this implies 

(■to 

p + 

and hence the claim. 



TT / i °( t o-s) 2 [c^l p+1 ](s,io-s)^<Flux(-l,t ) <E , 



□ 



IV. BOUNDEDNESS 



Using the light-cone estimate established in Proposition [T] the boundedness of ijj in the region 
K = K~ n{— 1 < t < 0} can be shown by a simple argument going back to Pecher Pec74 ] . 

Theorem 1. Let ip € C 2 (K) be a solution of (|13p with 2 < p < 5 and initial data supported on J 
satisfying the estimate (|14p . If c > is uniformly bounded and dtc < on K then ip is uniformly 
bounded. 

Proof. Let ipo G C 2 (K) be a classical solution of the homogeneous equation Oipo = with the 
same initial data as rp. Then, for a fixed (t, r) G -RT, it holds that 

1 rt—s+r 

\ip - ipo\(t,r) < — / / yc(s,y)|'0| p (s,2/)dyds. 
2r y_i 7|t- s _ r | 

Due to the fact that 2 < p < 5 there exists a q with 3/2 < g < (p + l)/(p — 1). Changing variables 
and applying Holder's inequality yields 

t rt—s+r 1 rt+r fb 



1 

2r 



1 J\t-s-r\ 



yc(s,y)\tp\ p (s,y)dyds 



1 

2r 



(to — s)c(s,to — s)\ip\ p (s , to — s)dsdto 



< 



2r 



t+r r b 



(to - s) 2 c q (s,t - sM pq (s,to - s)dsdt 



t+r rb 



9-2 



(to — s)i- 1 dsdto 
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(15) 
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where the abbreviations a := max{t— r, r— t— 2} and b := [iQ+(£— r)l /2 were introduced. Consider 
the first integral. Since q > 1 and c is bounded, so is c 9_1 . Furthermore, < pq — (p + 1) < q, so 
that 



(to - s)V(s,to - s)|VPMo - s)dsdt 

l 



/ \ 11 rt+r i-b 

<C[ sup ||^(r,-)IU«) / / (t -s) 2 c(s,t -s)\-i(j\ p+1 (s,t -s)dsdt , 

\-l<T<t J Ja J-l 

where < 7 < 1 is such that 7a = pq — (p + 1). The integration variable to takes values to > a 
which ensures — 1 < b < to- Thus, the integral in s can be estimated by virtue of Proposition CD to 
yield 

rt+r rb rt+r 

/ / (t - s) 2 c(s,t - s)\ip\ p+1 (s,t - s)dsdt < / (p + l)E dt < 2rCE . 

Ja J — 1 Jt—r 

The second integral in ([151) can be calculated directly to give 

f t+T f b q-2 n-l f t+V 2,-3 

/ (t - s) *-i dsdt < ^— - / (1 + to) ^ dt < 2rC 

J Ob J — 1 j. <J t — V 

because q > 3/2. To sum up, the estimate 

\^-H(f,r)<§-( sup U(T,-)\\ L ~Y(2rEo)*(2r)^ =CE*( sup \\^(r, -)\\ L ^) 
*r \-i< T <t J \-i<T<t / 

holds true for any (t, r) € if. But since the solution ipo of the homogeneous equation is clearly 
bounded and < 7 < 1 this estimate implies the boundedness of ijj itself uniformly on K, for the 
function 1 1— ► sup_ 1<r<i \\ip(T, •)||l° c is continuous and finite at t = —1. □ 

With the function ^ bounded on K the decay of <ft in the future of I follows immediately. 

Corollary 1. Let <fi be a classical solution of the wave equation (jlip for 3 < p < 5 with initial 
data 4>o G C 2 (R+) and </>i G C 1 (M+) given at t = 1 and which exists globally towards the future. 
Assume that the support of <po and <p\ is contained in the compact interval I C [0, a p [. Then there 
is a constant C G such that <p satisfies the decay estimate 



\<f>(t,r)\ <t; ttt ^ (16) 

Irv yi ~ (l + t + r)(l + t-r)P- 2 v 1 



for all t > 1 and < r < t. 



Proof. Given such a solution of class C 2 , it was shown in section [II] that its conformal transfor- 
mation ip = $>z(rL~ 1 (f)) is a classical solution of the wave equation (|13p on the future of the curve 
H in K~ , cf. Figure [U with 

" = — - — ^(nuvf- 1 



(p-2) 2 

according to equation ([T2]) . Now, for u < n < 

n — v 



$>*(riuv) 



1 1 1 
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so that c is bounded on the future of H in K owing to p > 3. Furthermore, 



1 $*(£luv 



,2 



p— 3 p— 3 



which implies <9(C < on IT - again by reason of p > 3. Moreover, as detailed in section IIIIl the 
assumptions on the support of (fto and 4>\ guarantee that J is compactly contained in {—1} x [0, 1[ 
and that the estimate (|14p holds. Thus, Theorem [1] applies and yields boundedness of tp on K. 
Since ip is also bounded on the compact region S it follows that if) is bounded on the whole future 
of H in K~, say \ip\ < C. But then 

m,r)\=n(t,rW^(t,r)\<Cn(t,r) 

on the whole future of I in K + , while 

1 



«(t,r) 



l-(^ 

u 



-2 ' 



Since in this region u > 1 and v > 1 — q p and outside, for v < 1 — q p , the solution vanishes 
identically, the bound for (j)(t,r) can be written in the regularized form (|16p . 

□ 



V. DISCUSSION 

As already noted in the introduction, the conformal transformations which we consider for the 
radial problem are not conformal for the full problem, beyond the spherical symmetry, unless p = 3. 
However, for the full problem with p > 3 we can still apply the conformal (p = 3) transformation 
and obtain a global pointwise decay l/(uv). The decay rate is then not optimal but can act as 
a prerequisite for a more refined asymptotic analysis. Indeed, that decay rate is sufficient for the 
solution to become small at large time s in th e sense which allows a perturbative analysis with 
methods similar to those developed in Szp08| . We want to address this issue in a forthcoming 



publication. 

It seems, at least as far as the spherical symmetry is concerned, that our method can be applied 
in higher than n = 3 odd dimensions. 

Also the linear wave equations with strong positive potentials 

d 2 t <\) -Acj) + V(x)<f> = 

having prescribed decay at spatial infinity V(x) ~ still lack a sharp pointwise decay estimate 

and seem to be treatable with our method after minor modifications. Such equations have a positive 
definite energy, too, and can be conformally transformed to a form analogous to (|13p with a regular 
positive function c. This idea shall be addressed in another publication. 
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